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ABSTRACT 



We discuss the nature of the brightest mid-IR point source (which we dub Object X) 
in the nearby galaxy M33. Although multi- wavelength data on this object have existed 
in the literature for some time, it has not previously been recognized as the most 
luminous mid-IR object in M33 because it is entirely unremarkable in both optical and 
near-IR light. In the Local Group Galaxies Survey, Object X is a faint red source visible 
in VRI and Ha but not U or B. It was easily seen at JHKg in the 2MASS survey. It 
is the brightest point source in all four Spitzer IRAC bands and is also visible in the 
MIPS 24 fim band. Its bolometric luminosity is ~ 5 x 10^ Lq. The source is optically 
variable on short time scales (tens of days) and is also slightly variable in the mid-IR, 
indicating that it is a star. Archival photographic plates (from 1949 and 1991) show no 
optical source, so the star has been obscured for at least half a century. Its properties 
are similar to those of the Galactic OH/IR star IRC+ 10420 which has a complex dusty 
circumstellar structure resulting from episodic low velocity mass ejections. We propose 
that Object X is a M > 3OM0 evolved star obscured in its own dust ejected during 
episodic mass loss events over at least ~half a century. It may emerge from its current 
ultra-short evolutionary phase as a hotter post-RSG star analogous to M33 Var A. 
The existence and rarity of such objects can be an important probe of a very brief yet 
eventful stellar evolutionary phase. 

Subject headings: stars: evolution, mass-loss, winds, outflows — galaxies: individual 
(M33) 



1. Introduction 

Explaining mass loss from massive stars, especially episodic mass loss in evolved massive 
stars, is one of the outstanding problems in stellar evolution theory. Systematic studies of rare, 
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luminous, d usty massive stars in nearby galaxies were carried out in parallel by iBonanos et al 



2009 . l20ld) for the La, r g;e an d Small Magellanic Clouds (LMC and SMC), and [Thompson et al 



20091 ) and iKhan et alJ (|2010|) for M33, NGC 300, M81, and NGC 6946. Luminous blue vari- 



ables (LBVs), supergiant B[e] (sgB[e]), some Wolf-Rayet stars (WRs) and red supergiants (RSGs) 
contribute to this class of objects, and their rarity implies they are a very short but perhaps crit- 
ical stage in the evolution of massive star s . Given that the LBV mass loss mechanism is poorly 
understood (jHumphreys &: Davidsonlll994l : ISmith Sz Qwockil 120061 ) and that th ere is mounting ev- 



idence for some core collapse supernovae with dust ens hrouded progenitors (iPrieto et al 



Thompson et al.l |2009| ) or very recent mass loss episodes (jGal-Yam et al.l 120071 : ISmith et al 



2008 



200e 



and references therein), a census of these stars in nearby galaxies is vital for understanding the 
mass loss mechanisms and final stages of evolution of massive stars. 



Bonanos et al.l (|2009l . l20ld ) cross-matched massive stars with known spectral types with the 
SAGE and SAGE-SMC photometric databases, which resulted in a multi-band photometric catalog 
from 0.3 to 24/xm providing spectral energy distributions (SEDs) for most classes of hot and cool 
stars. These works showed that in the LMC and the SMC, the LBVs, RSGs, and sgB[e] stars are 
some of the most-luminous mid-IR sources, due to the combination of high actual luminosities with 
surrounding dust from recent episodic mass ejections. However, these studies required the presence 
of a bright optical, spectroscopicall y classified source. The opp osite approa c h, foc using on bright 
mid-IR sources, was undertaken by [Thompson et al.l (|2009l ) and lKhan et al.l (j2010l ). 



The [Thompson et al.l (120091 ) re-analysis of the iMcQuinn et al.l (j2007l ) mid-IR images of M33 
focused on detecting extremely red (optically thick even at 3.6//m) analogs of the progenitors of 
SN 2008S and the NGC 300 optical transient, finding that such analogs are extremely rare. Less red, 
but more luminous, stars are still rare, and o nly in some cases d o they correspond to LBVs and other 
spectroscopically classified evolved stars. In lKhan et al.l (j2010l ). we carried out a systematic mid-IR 
photometric search for self-obscured stars in four galaxies: M3 3 , NG C 300, M81, and NGC 6946. 
In particular, we confirmed the conclusion of [Thompson et al.l (|2009l ) that stars analogous to the 
progenitors of SN 2008S and NGC 300 transients are truly rare — there may be as few as ~ 1 per 
galaxy at any given moment. This result empirically supports the idea that the dust-enshrouded 
phase is a very short-lived phenomenon in the lives of some massive stars. 

At this point, we set out to examine the mo st l uminous mid - IR so urces in M33 that had not 
been considered by either [Thompson et al.l (|2009l ) or lKhan et al.l (|2010l ). In this paper, we discuss 
the nature of the brightest mid-IR star in M33. We started by examining the brightest objects 



off th e primary stellar locus in the mid-IR CMD of M33 (Figure [H adopted from [Thompson et al 



20091 . also see Figure [2]). The brightest source is a known compact, young star cluster IC133. It 
is spatially unresolved by Spitzer, but immediately recognizable as a star cluster. It is both too 
luminous to be a single star [L^ > 10 ''Lq) and its mid-IR SED peaks at a temperature (~ 30K) 
that is too cold for dust associated with material ejected by a single star (because the required mass 
oc L^T^^). The second brightest source, the brightest mid-IR star in all of M33, is the subject 
of this paper. We caU this star, located at R.A.=1^33™24n and Dec.=+30°25'34'.'8 (J2000.0), 
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"Object X". 

Figure [3] shows its location in M33, and Figure H] shows the F-band through 24 /um band 
images of its surrounding region based on the data we describe in Section [2j Section [3] presents the 
SED and simple models of it. Section [H discusses the possible explanations of the physical nature 
of Object X and its implications for understanding the late-stage evolution of the most massive 
stars. 



2. Data 



The optica l ythotometTy (UBy RI and Ha) was measured from the Local G roup Galaxies 



Surve y images (jMassev et al.l 120061 ). The catalog of M33 point sources published bv iMassev et al 



(|2006l ) does not include Object X as a point source because it did not s atisfy their cr iterion that 
it is >4 Sigma above the background. We used DAOPHOT/ALLSTAR (|stetsonlll99a ) to identify 
point sources in the UBVRI bands. Object X is identified as a relatively faint and red point source 
in the VRI bands with a > 3cr detection. However, it is not detected in the U and B bands. We 
estimated 3a upper limits on its luminosity in these bands using the APPHOT/PHOT package. 
The measured magnitudes and the limits were transfor med to the V e ga-ca librated system using 
zero-point offsets determined from the bright stars in the lMassev et al.l (|2006l ) catalog of M33 point 



sources. 



The near-IR (JHKs) images were taken from the 2MASS survey (ISkrutskie et al.ll2006l ) an d the 



calibrated JHKs magnitudes were obtained from the 2MASS All-Sky Catalog of Point Sources (ICutri et al 
20o3). 



For the mid-IR photometry, we used the six co-add ed e pochs of data from iMcQuinn et al 



(120071) as processed and used bv [Thompson et al.1 ((20091) and iKhan et al.l tom ). For the MIPS 
24 fiin band, we downloaded the Post-Basic Calibrated Data (PBCD) from the Spitzer archive 
(Program 5, PI: Gehrz). Although we examined the MIPS 70 /im and 160 /im images for the 
target region as well, we were unable to obtain reliable flux measurements due to the nearby 
H II region being extre mely luininou s in these two bands. We measured the magnitudes using 
DAOPHOT/ALLSTAR (|stetsoJl99a ^. The PSF-magnitudes obtained with ALLSTAR were trans- 
formed to Vega-calibrated magnitudes using aperture corrections derived from bright stars using 
the APPHOT/PHOT package. Figures E] and [2] show the mid-IR CMDs of M33. 

The measured fluxes are reported in Table 1, and IZ-band through 24 ^m band images of 
Object X are shown in Figured! 

We also examined the photographic plates from the 48 inch Oschin Schmidt Telescope used 
in the STScI Digitized Sky Survey. The first image (epoch 1949) was a POSS-I Red Plate, while 
the second (epoch 1991) was a POSS-II RG610 filter. Both of these are roughly comparable to the 
modern standard -R-band. Neither image shows an optical source at the location of Object X. These 



- 4 - 



images are much shallower than the Local Group Survey (jMassev et al.ll2006l ). Based on the faintest 



-Riband USNO stars in this region and calibrating them by their counterparts in iMassev et al 



(|2006l ). we estimate an upper limit of R> 20.0 magnitude for Object X in 1949. This is consistent 
with our detection of Object X at R~ 21.6 (epoch 2001). The epoch 1991 image also does not show 
a source, but we were unable to converge on a self-consistent calibration for the region, and so will 
not discuss this image further. Figure [5] shows these historical images of this region. 

Object X is identified as an optically variable point source in the Canada- France-Hawaii Tele- 
scope (CFHT) photometric survey of M33 (jHartman et al.ll2006l ). Figure [6] shows the r' and i' band 
lightcurves of Object X from this survey. The correlated short-term variability of ~ 0.4 magnitude 
(fractional variability of '^45%), definitively indicates that it is a single stellar object rather than 
multiple objects blended together. If we were to dilute the source with additional non- variable 
sources of comparable total luminosity, the observed ~ 0.4 magnitude variability would be unphys- 
ical. This is significant given the lack of archival HST images for this location. 



McQuinn et al.l (j2007l ) collected six epochs over two years and identified Object X as a vari- 



able source in the mid -IR, although for this paper, we use the mid-IR lightcurves generated by 
Thompson et al.l ((20091). Figure [7] shows the 3.6 /im and 4.5 ^m band lightcurves. Object X shows 
correlated variability of about ~0.15 magnitude in both bands (fractional variability of ~15%). 
This suggests that the dust opacity and geometry surrounding Object X are not undergoing any 
rapid evolution. 



3. Analysis 



Figure [8] shows the spectral energy distribution (SED) of Object X. The low optical luminosity 
means that Object X is unlik ely to be a YSO emb edded in a dusty disk, and if we model the 
SED using the YSO models of Robitaille et al. (2006), we find no good, or even plausible fit^. In 
essence, <1% of the luminosity is emitted in the optical, while even edge-on YSO disk models tend 
to scatter more optical light into the observer's line of sight. 

The alternate is that we are examining a self-obscured star. For this case, we fit th e SED of 
the source using DUSTY dlvezic Elitzurlll997l : [lvezic et al.lll999l : lElitzur fc Ivezicl[200ll ') to model 



the radiation transfer through a spherical dusty me dium surroundi i ig a s tar. We considered models 
using either the graphitic or silicate dust models of lDraine fc Led (j 19841 ). We distributed the dust 
either as a shell, with ratios of outer to inner radii of either 2:1 or 4:1, or with the standard DUSTY 
wind model. The density distributions inside the shells were assumed to be p oc 1/r^. The models 
are defined by the stellar luminosity, L^,, stellar temperature, T^,, the V-band optical depth ry, and 
the dust temperature at the inner edge of the dust distribution, T^. The stellar luminosity and the 



^We fit the SED using tiie online tool l|Robitaille et al.lbOOTI ') at 
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dust temperature together determine the radius of the inner edge of the dust distribution. To fit 
the data we tabulated the DUSTY models on a grid of stellar temperatures, optical depths and 
dust temperatures. 

We find that two dust types that can fit the SED equally well but with very different preferred 
stellar temperatures. The graphitic models prefer cool stars, ~ 5000 K, while the silicate 
models prefer hot stars, ~ 20000 K. The stellar luminosity is L^, ~ IO^'^Lq to lO^'^L© with 
the silicate models favoring the somewhat higher luminosities (see Figure [9]). The three different 
dust geometries produce similarly good statistical fits. For the fits we assumed minimum flux 
uncertainties of 20% and included the optical upper limits as constraints, leading to — 43 for 
N(i(,f = 9 for the best fits. Each model was a modestly poor match to some part of the SED, but 
this is not very surprising given our use of a simple spherical model and a discrete sampling of its 
parameter space. A clear detection (or not) of the silicate dust feature near 10 //m would determine 
the dust composition. The typical optical depths and inner/outer dust temperatures were Ty — 8.5 
and Td ~ 500/150 K for the graphitic models, and Ty — 11-5 and To ~ 1200/300 K for the silicate 
models. Much of this difference is created by the higher scattering opacities of the silicate dust. 
The inner edge of the dust lies at ~10^^'^ cm for the graphitic models and ~10^^'^ cm from 
the silicate models. 



Based on the CO and H I data from ([Gratier et al.ll2010l ). the estimated column density at 
the location of Object X is ~ 7 x 10^^ atoms/cm^ (P. Gratier, private communication), where the 
higher number comes from using the CO to estimate the amou nt of H?. For a st andard Galactic 



dust to gas ratio of E{B — V) = 1.7 x 10 mag cm /atom (iBohlin et al.lll978l ). this implies a 



maximum foreground extinction in M33 of roughly E{B — y) ~ 1.2 magnitude. Adding this to the 
Galactic extinction changes the DUSTY models very little, modestly reducing the optical depths 
(Ar ~ 3 — 4) and slightly raising the stellar luminosity (by ~20%). 



If we search the Padova stellar models (iMarigo et al.l 120081 ) for stars with such luminosities, 
there are examples with both the high and low stellar temperatures, although the higher tempera- 
ture stars of the silicate models are preferred. The stellar masses are always > 3OM0, but without 
an independent constraint on the stellar temperature it is difficult to say more. In many, but not 
all cases, the stars have lost significant fractions of their initial masses. 

We also note that while stars ordinarily do not have significant Ha emission. Object X is a 
strong Ha source. We determine this by using difference imaging methods to match the Ha image 
to the flux scale and PSF structure of the i?-band image and then subtract. All the normal stars 
disappear to leave us only with the st ars having significant Ha emission. I n practice, we use the ISIS 
ir nage subtrac t ion so ftware package (lAlard h LuptonI Il998l : lAlardI 12000) following the procedures 
of iKhan et al.l (j2010l ). Figure [TOl shows the result of this procedure for the region around Object X. 
It is clear from the subtracted image that Object X is a strong Ha emission source. 
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Discussion 



An evolved star can cloak itself in dust through two broad mechanisms, winds and mass ejec- 
tion. The classic examples of wind obscured stars are the AGB stars. In the most extreme cases, 
such as the progenitors of SN 2008S and the 20 08 NGC 300 transient, the wind can be opti- 
cally thick even in the mid-IR (jPrieto et al.l 120081). Shells formed by impulsive mass ejections are 



seen around many evolved massive stars (e.g . . iHurtiphrevs &: DavidsonI 1994 : [Humphreys et al 



1997 



Smith Sz Frew 



201C 



Gvaramadze et al.l l20lJ[Wachter_et_al 



2010|) 



and they are distin- 
guished by the frequency of the ejections, their mass, and their velocities. The most famous 
example is the "Great Eruption " of r? Carina in the 19th cent ury, which ejected ~ IOMq of ma- 



terial at velocity ~ 600 kms ^ ( Humphrevs Sz DavidsonI Il994l ). At the other end of the velocity 
spectrum ar e the QH/IR stars such as IRC+10420, where the ejection velocities are closer to 
~ 50 kms-i (|Tiffanv et al.ll201ol ). 



Our best constraint on these possibilities comes from the time variability of the source. The 
absence of the source in 1949 means that the star has been obscured for at least 60 years, but the 
SED models require fairly hot dust close to the star. If we assume that the material obscuring 
the star today is the same as that obscuring the star in 1949, we get characteristic velocities of 
the order of ~ 114 kms~^ and ~ 36 kms~^ for the graphitic and silicate models, respectively. The 
outer edge then requires a velocity twice as large for a shell with a 2:1 ratio between the inner and 
outer edges. These velocity scales favor the slower ejections of the OH/IR stars over the higher 
velocities of the giant LBV eruptions like rj Carina. The significant optical and mid-IR variability 
then favors a patchy evolving ejecta over a steady wind. The extended period of the obscuration 
comes close to requiring multiple ejections, approaching the limit of an unsteady wind. It would 
be interesting to determine if the source was visible in still older archival plates, if any exist, since 
it would be a relatively bright R ~ 16 mag source without the dust. 

Overall, Object X seems to most closely resemble the cool, hypergiant, stars such as IRC+10420 
and Var A. The three sources have comparable bolometric luminosities of ~ 10^' ''Lq, although their 
SEDs, shown in Figure [HI differ because Object X is more obscured. In the case of IRC+10420, the 
star is surrounde d by a series of shells and seems to have undergone a period of high mass loss for the 
last ~ 600 v ears (IHumphrevs et al .11 199 71 ) that ended 50-100 years ago leading to its current optical 
brightening ( Blocker et al. 19991 ). Var A, on the other hand, seems to have had a brief period of high 
mass loss rates over the last ~ 50 years and is now re-emerging in the optical (IHumphreys et al 



20061 ). Object X would seem to be intermediate, requiring a longer period of heavy mass loss than 
Var A, but perhaps less than IRC+10420. If Object X is a true analog to IRC+10420 or Var A, 
then it may conceivably emerge from its current self-obscured state over the next few decades since 
the optical depth of an expanding shell drops as r oc We cannot make this statement with any 
degree of certainty, as each of these objects are going through tumultuous evolutionary phases that 
are both unique and poorly understood. Follow-up observations of this source at all wavelengths 
should begin to reveal its long term evolution and may provide us with unexpected surprises. Deep 
optical and IR spectra of Object X will enable to us to determine the spectral type of the central 
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star and place it on an H-R diagram. 



Stars like Object X are extremely rare. As a unique star in a brief yet eventful evolutionary 
state, we encourage further study of this interesting object. Furthermore, there is mounting ev- 
idence that mass l oss from massive stars may be dominated by impulsive transients rather than 



steady winds (e.g., ISmith &: Owocki 



200 



pa rticularly with the downward revision of mass-loss 



Fullerton et al. 



20061 ). If so, the period in which mass loss is most impor- 



rates in recent years (e.g., 
tant will als o tend to be t h e peri od w hen the star is most o b scure d. Surveys of massive stars such 
as those by iMassev et all (|2006l ) and iBonanos et al.l ((20091, I2OIOI ) , which focus on bright optical 
sources, will miss Object X and similar sources that may best probe the frequency and duration of 
these mass loss episodes as well as the amount of ejected mass. Characterizing these stars and the 
amount of mass loss by massive stars clearly requires systematic surveys in the mid-IR as well as 
the optical. 
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0.049 



0.96 Mpc (IBonanos et al.l |2006| ) and corrected for Galactic extinction o f E{B 
dSchlegel et al.lll998l). The UBVRI images are from iMassev et al.l (l20odl. the JHK , photome- 
try is from ICutri et al.l (j2003l ) , the IRAC bands photometry is from [Thompson et al.l (120091 ) , and 
the MIPS 24 //m band image from the Spitzer archive (Program 5, PLGehrz). 



-9- 



Bonanos, A. Z. et al. 2009, AJ, 138, 1003 

— . 2010, AJ, 140, 416 

Cutri, R. M. et al. 2003, 2MASS All Sky Catalog of point sources., ed. R. M. Cutri et al. 
Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89 
Elitzur, M. & Ivezic, Z. 2001, MNRAS, 327, 403 

FuUerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 
Gal- Yam, A. et al. 2007, ApJ, 656, 372 
Gratier, P. et al. 2010, A&A, 522, A3+ 

Gvaramadzc, V. V., Kniazcv, A. Y., & Fabrika, S. 2010, MNRAS, 405, 1047 

Hartman, J. D., Bcrsicr, D., Stanek, K. Z., Beaulieu, J., Kaluzny, J., Marquette, J., Stetson, P. B., 
k Schwarzenberg-Czerny, A. 2006, MNRAS, 371, 1405 

Humphreys, R. M. & Davidson, K. 1994, PASP, 106, 1025 

Humphreys, R. M. et al. 1997, AJ, 114, 2778 

— . 2006, AJ, 131, 2105 

Ivezic, Z. & Elitzur, M. 1997, MNRAS, 287, 799 

Ivezic, Z., Nenkova, M., & Elitzur, M. 1999, arXiv:astro-ph/9910475 

Jones, T. J. ct al. 1993, ApJ, 411, 323 

Khan, R., Stanek, K. Z., Prieto, J. L., Kochanek, C. S., Thompson, T. A., k Beacom, J. F. 2010, 
ApJ, 715, 1094 

Marigo, P., Girardi, L., Bressan, A., Groenewegen, M. A. T., Silva, L., & Granato, G. L. 2008, 
A&A, 482, 883 

Massey, P., Olsen, K. A. G., Hodge, P. W., Strong, S. B., Jacoby, G. H., Schlingman, W., & Smith, 
R. C. 2006, AJ, 131, 2478 

McQuinn, K. B. W. et al. 2007, ApJ, 664, 850 

Prieto, J. L. et al. 2008, ApJ Letters, 681, L9 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJ Supplement, 169, 328 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJ Supple- 
ment, 167, 256 



-10- 



Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Skrutskie, M. F. et al. 2006, A J, 131, 1163 

Smith, N., Chornock, R., Li, W., Ganeshalingam, M., Silverman, J. M., Foley, R. J., Filippenko, 
A. v., & Barth, A. J. 2008, ApJ, 686, 467 

Smith, N. k Frew, D. J. 2010, arXiv:astro-ph/1010.3719 

Smith, N. & Owocki, S. R 2006, ApJ, 645, L45 

Stetson, R B. 1992, 25, 297 

Thompson, T. A., Prieto, J. L., Stanek, K. Z., Kistler, M. D., Beacom, J. F., k Kochanek, C. S. 
2009, ApJ, 705, 1364 

Tiffany, C, Humphreys, R. M., Jones, T. J., & Davidson, K. 2010, AJ, 140, 339 

Wachter, S., Mauerhan, J. C, Van Dyk, S. D., Hoard, D. W., Kafka, S., & Morris, R. W. 2010, AJ, 
139, 2330 



This preprint was prepared with the AAS LM^^X macros v5.2. 



- 11 - 



14 



1 — I — I — I — I — I — I — I — I — I — I — r 



T 



IC133 » 
Object X» 

IRC+10420 



12 



10 



o Var A 



-8 



I ' I ' ■ j I ' I ■ .1 . 



J I I I I I I I I I I L 



-0.5 0.5 1 

[3.6]-[4.5] 



1.5 



Fig. 1. — The 4.5 Atm a bsolut e magnitude vs. the [3.6] — [4.5] color mid-IR CMD of M33 adopted 
from iThompson et al.l ()2009l ). Object X (red circle), the star cluster IC133 (cross and circle) 
and Var A (open circle around dots) are marked. For comparison, we a lso show the position of 
IRC+10420 (blue starred symbol) for a 5 kpc distance (| Jones et al.lll993l ). The extremely bright 
stars near color ~ are the foreground stars and the stars on the red branch at M4.5 ~ —10.5 are 
M33 AGB stars. There are no other stars in M33 that are even remotely similar to Object X in 
these bands. 
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Fig. 2. — Same as Figured] but for 8.0 /xm absolute magnitude vs. the [4.5] — [8.0] color. Object X 
stands out in this combination of mid-IR bands as well. 
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Fig. 
Dec. 



3.~ The IRAC 3.6 fim image of M33 showing the location of Object X (R.A. = 1^33'"24n, 
= +30°25'34'.'8; J2000.0) with an arrow. 
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Fig. 4. — Multi-band images of the region around Object X (marked b y the circles). The optical 
images were taken from the Local Group Survey (jMassey et al.l 120061 ) . the near-IR images were 
taken froi n 2MASS (ISkrut skie et al.ll2006l ). the mi d-IR images are from six co- added epochs of the 
data from iMcQuinn et al. (2007,) as processed by [Thompson et al.l (j2009l ). and the MIPS 24 ij,m 
band image is from the Spitzer archive. 
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Fig. 5. — Historical images of the region near the location of Object X (marked by the circles) over 
the last ~ 60 years implying that Object X has remained obscured at least over this period of time. 
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Fig. 6. — Optical variability of Object X in the r' and i' bands from iHartman et al.l ()2006l ). 
Object X shows correlated variability of about ~0.4 magnitude (fractional variability of ~45%) in 
both bands. 
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Fig. 7. — Mid-IR variability of Object X in the 3.6 fim and 4.5 //m IRAC bands. It shows correlated 
variability of about ~0.15 magnitude (fractional variability of ~15%) in both bands. 
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Fig, 



The SEP of Object X as compared to tha t of IRC+10420 (jJones et all Il993l : 



Humphreys et al.l 119971 ) and Var A ([Humphreys et al.l 120061 ) . Both comparison SEDs haye been 



constructed using data from multiple epochs (IRC+10420: UB from 1972, VRI from 1992, the rest 
from 1996; Var A: optical from 2000-01, near-IR from 1997, the rest from 2004-05) and these sources 
are known to be v ariable. We also note the large amount of extinction toward IRC+10420 {Ay — 6, 



Jones et alJll993l ). correcting for this extinction would make it more similar to Var A than Ob- 



ject X. However, it is uncertain how much of the extincting material is associated with that star 
rather than simply being along the line of sight, and therefore we do not correct for the extinction 
here. 
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Fig. 9.— The results of DUSTY fits to Object X SED (black points). The heavier sohd (dashed) 
lines show the model DUSTY SEDs for the source for graphitic (silicate) dusts. The lighter solid 
(dashed) lines show the black-body stellar SEDs. The graphitic model has T* = 5000 K and 
= 10^'^Lq with Ty = 8.5 and an inner edge dust temperature of 500 K corresponding to an 
inner shell radius of 2.6 x 10 cm. The silicate model has = 20000 K and = IO^-^Lq with 
Ty = 11.5 and an inner edge dust temperature of 1200 K corresponding to an inner shell radius of 
6.4 X 10^^ cm. These models have a 2:1 ratio between their inner and outer radii, but the 4:1 and 
wind models (not shown) look very similar. 
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Fig. 10. — The Ha, R-band, and Ha— R subtracted images of the region near Object X. In the 
subtracted image, almost all the objects vanish, leaving behind the nearby H II region, Object X, 
and a small number of additional, mostly extended, Ha sources. 



